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PREFACE 


This  report  describes  work  done  during  1984  by  Dr  J.  F.  Mahoney, 
investigator,  from  the  Department  of  Industrial  and  Systems  Engineering,  the 
University  of  Florida,  Gainesville,  Florida  32611,  under  Contract  No.  F08635- 
83-C-02n2  sponsored  by  the  Air  Force  Armament  Laboratory  (AFATL),  Armament 
Division,  Eglin  Air  Force  Base,  Florida  32542.  The  program  manager  was 
Mr  Daniol  A.  Me  Inn  is  (DLYW). 

The  work  provides  the  mathematical  and  statistical  techniques  of  the 
methodology  developed  in  the  four  series  of  reports  entitled  "Variability  of 
Measures  of  Weapons  Effectiveness,"  Volumes  I,  II,  III,  and  IV,  in  the  form 
of  equations  that  can  be  coded  for  high  speed  computers. 
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SECTION  I 


INTRODUCTION 

This  report  summarizes  the  important  equations  and  calculations  which 
were  developed  in  the  four-volume  research  report  "Variability  of  Measures  of 

* 

Weapons  Effectiveness"  written  primarily  by  B.  S.  Sivazlian,  [21, [3], [4],  and 
[5].  Section  II,  III,  IV,  and  V  correspond  on  a  one-to-one  basis,  to  Volumes 
I,  II,  III,  and  IV  of  the  above  mentioned  report.  The  objective  is  to  present 
equations  for  the  estimation  of  probability  of  kill  and  its  variance  for 
targets  which  are  only  sensitive  to  blast  or  to  fragments  when  there  is 
delivery  error  and  when  there  is  not. 
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SECTION  II 


FRAGMENTS  SENSITIVE  TARGETS  IN  THE 
ABSENCE  OF  DELIVERY  ERROR 


1.  Objective 


The  equations  necessary  to  compute  E[P)(f]  and  Var[P|^f:]  in  the  absence  of 
delivery  error  are  presented.  Two  approaches  are  used. 


2 .  Input  Parameters 

Numerical  values  of  the  following  must  be  provided: 


R  =  minimum  estimate  of  the  weapon  radius  in 
the  direction  of  range,  [ft]. 

R^  =  maximum  estimate  of  the  weapon  radius  in 
^2  direction  of  range,  [ft]. 

R  =  minimum  estimate  of  the  weapon  radius  in 
the  direction  of  deflection,  [ft]. 

R  =  maximum  estimate  of  the  weapon  radius  in 
the  direction  of  deflection,  [ft]. 

X  =  impact  coordinate  in  the  direction  of 
range  where  the  target  is  at  x=0,  [ft]. 

y  =  impact  coordinate  in  the  direction  of 
deflection  where  the  target  is  at 
y=u,  [ft]. 


3 .  Spe£i al  Functions 

One  must  have  access  to  P(XlY)  which  is  the  chi-squared  probability 
function.  In  SAS  this  is  supplied  by  PROBCHI(X,Y) .  The  related  function 
Y)  =  1  -  R(X!y)  may  be  found  in  Table  26.7  of  Abramowitz  and  Stegun  [1]. 
for  cn*'r.kinij  purposes  P(7|10)  =  .27456. 
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4.  Defined  Function 

The  following  function  is  now  easily  computed: 


_  [P(2-2|l)  -  P(2-2|l)]  B  exp{-  -o]  -  A  exp(-  -j) 
J(A.B.k) - +  - — - 


This  equation  is  a  recasting  of  (7)  of  [2],  It  was  obtained  from  (7)  by  first 
dividing  through  by  B-A  and  then  using  the  substitution 


*(x)  = 


The  prime  benefit  is  that  the  first  two  equations  of  the  next  section  are 
simpler  in  appearance.  These  two  equations  are  exactly  (8)  and  (10)  of  [2]. 


5.  Subjective  Estimation  Method 


ECP,,f]  =  J(R^  .  .  y)  0(R^  .  R,  .  x) 

2 


E[pL]  =  J(R^  ,  Ry  ,  y/2)  J(R  ,  R  ,  x/2) 
-^1-^2  1  2 


VarCP,,]  =  E[P^,]  -  (E[P^f])2 


3 


r. - 


4 


SECTION  III 


BLAST  SENSITIVE  TARGETS  IN  THE 
ABSENCE  OF  DELIVERY  ERROR 

1.  Obj ecti ve 

The  equations  necessary  to  compute  E[P|(f]  and  Var[P|(f]  for  blast 
sensitive  targets  in  the  absence  of  delivery  error  are  presented. 

2.  Input  Parameters 

Numerical  values  of  the  following  must  be  provided; 

Aj  =  minimum  estimate  of  the  radius  for  which 
kill  is  total ,  [ftl. 

A2  =  maximum  estimate  of  the  radius  for  which 
kill  is  total ,  [ft], 

Bj  =  minimum  estimate  of  the  radius  for  which 
there  is  no  kill,  fftl. 

B2  =  maximum  estimate  of  the  radius  for  which 
there  is  no  kill,  [ft]. 

These  numbers  conform  to 

0  «;  A^  <  A^  <  Bj  <  B^  “ 

3.  Annular  Regions  and  Special  Conditions 

Five  annualar  regions  are  identified  which  are  centered  about  the  point 
of  impact.  R  is  the  distance  in  feet  from  impact.  The  regions  are: 


Region 

a: 

0  <  R  <  A^ 

Region 

b; 

A^  <  R  <  A^ 

Reg  ion 

c : 

Ap  <  R  <  Bj^ 

Region 

d: 

Bj  <  R  < 

Region 

e: 

Bj,  <  R  <  ® 
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In  addition,  the  relations  among  Aj ,  A^,  ,  and  admit  to  four 

conditions.  They  are 


^2 »  ^1  ~  ~  ^ 


Condition 

1: 

=  A 

Condition 

2: 

<  A 

Condition 

3: 

=  A 

Condition 

4; 

<  A 

4.  Subjective  Estimation  Method 

2 

The  equations  for  calculating  E[P|(f],  E[P|^^],  and  Var[P|^f]  have  been 
extracted  from  [31.  The  table  given  below  acts  as  a  directory.  One  must 
first  choose  which  region  (a,b,c,d,  or  e)  is  of  interest  and  then  which 
condition  (1,?,3,  or  4)  applies.  Then,  from  the  table,  find  the  reference 
number  which  directs  one  to  the  appropriate  set  of  equations  which  are 
presented  after  the  table. 


Condition  1  (1) 


Condition  2  (4)  (5)  (6) 


Condition  3  (8) 


(9)  (10)  (11) 


Condition  4  (12)  |  (13)  \  (14)  (15)  (16) 


ECP,,]  -1 
-  1 

VarCP^j,]  =0 

EtPkbT  ■  (M)' 

Var[Pk^]  '  0 


Et^kbl  ■  “ 


«^kbl  =  » 


Var[Pkbl  ■  0 


«\bl  =  ' 


E[P',1  =  I 


VarrPki,]  •  0 


B-R 


A^-Ai 


1  n 


B-A,  ^  A^_ 


R-R 


A^“i 


(B-R)(R-A,)  A^- 


(B-A^)(A2-AJ  ^  A^- 


''P'T%bl  ■  Efkbl  -  ('^t'’kb’) 


■  (B-A^)  (B-A^) 


(7) 


(B) 


(10) 


VarCP^^]  ■  ErP^^l  -  (Er.Pk(,1)^ 

ECEkbl  ■  0 

ECE^bl  ■  0 

Var[P^^]  =  0 

E^kbl  '  1 

Eff^l  =  ' 

VarCPkt,]  =  0 

R  A  B«-A 

''■"■E%bl  =  ff’^kb''  -  (EEPk^l)^ 

R-A  Bj)-A 
^^^kb"'  "  B^-B^  ■  B^-B”  ~R-A 

,  B,-R  (R-A)(R„-R)  P  . 

^-'■^bJ  "  B,-R^  (Bp-B^)(B2-A)  ^B^-B,^ 

VarrPk.l  =  ErpJ,l  -  (EEPk,!)^ 

Ef-kbl  ■  " 


e[pL]  ■  0 


Var[P^^l  .  0 


(12)  E[P|^„]  .  1 

ECP^^]  ■  1 
Var[P^^]  =  0 


il) 


9  . 


(13)  ^  ^  (A2-A^)(B2-B^)  b'^-R 


-  (B^-R)^  fjTT  -  B^l 


2  , 


'  - 


(R-AJ 


2 


1 '  ®2”^1 

1 - — L 


(A2-A^)(B2-B^)  B^-A^ 


Va^CPkb^  =  -  (ECP,,])' 


(\) 


0^)  •  J  •  (aP^TTvV  '"  '" 


-  O-*!''  "’^*  '"4^’ 


2  ®2"^2i 


(R-A  )^  Bp-A^  (R-A^)^  B^ 

Et^kbl  ■  >  -  Tv^TvSTT  ’"  li;^;-  *  (A;,-A;)(b^--b7  ’"  it; 


VarrP.J  -  ECP^J  -  (E[P,,])^ 


^  ■  .V 


V  ^  -'1  -.■  -.«  VFM  —*  •  1.  -  V-'.  *\  ••.  'T. 


’  ^  ^ 


•  •  •  • 

•N 


(B^-R) 


i\) 


7.  ,  ®2'^2 


.2  ,_  ®2'^1 


(IB)  F[P^^1  =  ^(B^-B^)  ^  (A^-A^)(B2-R|  R_A^  ‘  R-A^ 


-  < V">  BpA^l 


B,-R 


B„-A„ 


R.-A, 


*  (vviB-iv  R?^ -  '" -TT^l 


»ar[P|,j,]  =  EfpJj  -  (E[P,J)'' 


(16)  EtP,,^]  =  0 


ECP^V  =  0 


kb- 


kb 


VarCPkb]  =  0 


k 
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SECTION  IV 


FRAGMENT  SENSITIVE  TARGETS  IN  THE 
PRESENCE  OF  OELIVERY  ERROR 


1.  Objective 

The  equations  necessary  to  computer  E[P|(f]  and  Var[P|(f]  for  fragment 
sensitive  targets  in  the  presence  of  delivery  error  are  presented.  Two 
approaches  are  used. 


?.  Input  Parameters 


Numerical  values  of  the  following  parameters  must  be  provided. 

R  =  minimum  estimate  of  weapons  radius  in  the 
direction  of  range,  [ft]. 

R  =  maximum  estimate  of  weapon  radius  in  the 
^2  direction  of  range,  [ft]. 

R  =  minimum  estimate  of  weapon  radius  in  the 
direction  of  deflection,  [ft]. 

R  =  maximum  estimate  of  weapon  radius  in  the 
^2  direction  of  deflection,  [ft]. 

a  =  minimum  value  of  the  standard  deviation  of  the  aiming 
error  in  the  direction  of  range,  [ft]. 

a  =  maximum  value  of  the  standard  deviation  of  the  aiming 
error  in  the  direction  of  range,  [ft]. 

o  =  minimum  value  of  the  standard  deviation  of  the  aiming 
■^1  error  in  the  direction  of  deflection,  [ft]. 

a  =  maximum  value  of  the  standard  deviation  of  the  aiming 
^2  error  in  the  direction  of  deflection,  [ft]. 

3.  Defined  Quantities 


The  eight  quantities  listed  below  must  be  calculated. 
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4 .  Subjective  Estimation  Method 
Calculate  the  following: 


ECl,] 


- i (X  .  X  )  -  a/  (Xj,,  -  X  ) 

/2  (R,  -  )(a^  -  a  ’'e  ^2  Xj  21  11 

^2  ^2  ^1 


2  “  °x7‘rwx-i) 


E[Iy] 


^  (Ry  -  Ry  )(%  -  %  )  ^2  ^  ^ 

^2  ^1  ^2  ^1 


+  -  Jn  - -22 - 1  j„  _Z  12  ] 

^  "yi  (1  ^  Y^j)  2  a  -  Y^^T 


ECP^f]  =  E[IJ  E[Iy] 


E[i2] 


2  "tr; 


Rx  MOv 
1  ^2 


1  j  2  p 

r%  ^22  .  -1  1  °*1  ^21  ,  -1  1 

”K,)  '  2  Hz'  2”  % 


2  p  2  . 

'^x  x2  ®x  x2 

^^2  ^12  -1  1  ^  *^1  ^11  .„-!  1  , 
- A - sin  —  + - ^ -  sin  y — I 

^  ^12 


v  J:-^  v:^:  :  >  -■  V  ^ 


• . '  - ' .  •', ', .  ■ . 


"«  *.  *.  V”  '.  -^.  \ 


/2 

%  )K 


<  -j 

r  ^2  '22  -1  _1 

[  ^  Sin  y 


O  y2 

1  ^1  ^21  .  -1  _1_ 
Y  -  P  sin  y 


^  o  ^  2  , 

^2  42  .  -1  J_  .  ^1  ^11  ^.„-l  _L.i 

r\  sin  V  ^  0  sin  y  J 

Yj2  ^ 


Var[P^^l  =  E[I^]  E[i2]  -  (£[2^^])^ 


2  2 

The  equations  for  E[I  1  and  E[I  1  that  are  given  above  are  mere 

X  y 

rearrangements  of  (23)  and  (24)  of  [4].  In  the  versions  just  presented, 
greater  use  of  the  X^-j  and  Y^j  has  been  made  and  also  identities  of  the  form 


-1 


/2 


X?.  -  1 


sin  y  ~ 


have  been  incorporated. 


5.  Taylor's  Series  Method 


Estimates  of  E[P|^fl  and  VarlPj^^l  may  be  computed  using  a  Taylor's  series 
approach.  First  calculate. 


R  .  K 

X2  +  Xj 


R  ^  R 
^2  " 
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a  .  o 

^2  VA 


^  =  2 


Var[R^] - ^ - 


Var[Ry]  = 


2 


VarCo^l 


VarCa^]  = 


2 


T? 


Then  compute 


and 


ECPkfl  = 


/  R^  +  2a^ 


VarCP^fl  = 


F-[P 


1  + 


2 


Var[R^] 


X 


Varfa^l 
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SECTION  V 


BLAST  SENSITIVE  TARGET  IN  THE  PRESENCE 
OF  DELIVERY  ERROR 


1.  Objective 

The  equations  necessary  to  compute  Var[P|^(jl  for  blast 

sensitive  targets  in  the  presence  of  delivery  error  are  presented. 

2.  Input  Parameters 

Numerical  values  of  the  following  must  be  provided. 


Aj  =  minimum  estimate  of  the  radius  for  which 
kill  due  to  blast  in  complete,  [ft]. 

A2  =  maximum  estimate  of  the  radius  for  which 
kill  due  to  blast  is  complete,  [ft]. 

=  minimum  estimate  of  the  radius  for  which  there 
is  no  kill  due  to  blast,  [ft]. 

=  maximum  estimate  of  the  radius  for  which  there 
is  no  kill  due  to  blast,  [ft], 

0  =  minimum  value  of  the  standard  deviation  of  the  aiming  error 

*1  in  the  direction  of  range,  [ft]. 

o  =  maximum  value  of  the  standard  deviation  of  the  aiming  error 
^2  in  the  direction  of  range,  [ft]. 

a  =  minimum  value  of  the  standard  deviation  of  the  aiming  error 
^1  in  the  direction  of  deflection,  [ft]. 

a  =  maximum  value  of  the  standard  deviation  of  the  aiming  error 
^2  in  the  direction  of  deflection,  [ft]. 
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4.  Special  Function 


One  must  have  access  to  P(XlY)  which  is  the  chi-squared  probability 
function.  In  SAS  this  is  supplied  by  PROBCHI(X,Y) .  The  related  function 
0(X|Y)  =  1  -  P(X|Y)  may  be  found  in  Table  26.7  of  Abramowitz  and  Stegun  [1]. 
For  checking  purposes  P(7|10)  =  .27456. 

5.  Additional  Oefined  Quantities 

The  quantities  X^(^),  X^(B),  X2(A),  X2(W) ,  Y^fA)  and  Y2(B’)  are  all  given 
by  infinite  series  which  are  of  the  form 

C[tj^  ^2  ^  tj^  +  ...  ] 

The  definitions  of  the  entering  items  in  this  series  differ  depending  on  which 
of  the  six  quantities  is  to  be  calculated.  The  series  all  converge  rapidly  so 
that  only  the  first  few  terms  of  the  series  need  be  generated. 

(a)  Calculation  of  X^CS). 


t^  =  P(cxl4n-2)  C 


n-1 


(2n-2)! 
(n-1) !  (n-i; 


n=l,2,... 


Xj(A)  =  C[tj 


where  N  is  the  number  of  terms  retained  in  the  series.  N  may  be  taken  as  the 
smallest  integer  for  which 
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tn  =  P(Bl4n.2)  C 


X^(B)  =  C[t^  *  ^ 


wh(>rt‘  N  is  the  si^allest  integer  for  which 


^  <  10-' 


(c)  Calculation  of 


C  = 


/  n 


4a 


3/2 


=  P(oi|4n-l)n 


n-1 


4n- 


(n-D!  (n-1)!  (2n-2)! 


X^(A)  =  C[t^  +  +...+  t, 


where  N  is  the  smallest  integer  for  which 


N  -ft 
^  <  10 

1 


'0)  Calculation  of  X5(B) 


tfi  "  P(6|4n-l)n  ^  (n-D!  (n-1  (2n-2) 

X^CR^  =  C[t^  +  +...+  tj^] 


where  N  is  the  smallest  integer  for  which 


^  <  10-6 
1 


(e)  Calculation  of  'f . 


C  = 


b/iT 


32  a 


5/2 


tn  *  P(“l4n+l)n  ^  n!  (n^l)i^(2n-l) !  ’  "  " 


Y^(A)  =  C[t^  +  tp  +...+  tj^' 


where  N  is  the  smallest  integer  for  which 


10 


-6 


(f)  Calculation  of  ^2^^) 


C  = 


32a 


b/fl” 


=  I rTCn-'llMLl)! 


=  l,^. 


Y^(R)  =  C[tj  +  t^  +...+  tj^ 
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where  N  is  the  smallest  integer  for  which 


10 


-6 


Items  (a)  through  (f)  given  above  are  rearrangements  of  (54)  and  (55)  of 
[5]  for  particular  parameter  values. 


6.  Taylor's  Series  Method 


E[P 


RXj(B)  -  AX^(A)  -  X5>(B)  +  X^(A) 

(¥-A)o^Oy 


Then  calculate  the  quantities 


R 


1 


B-A  a^ay(B-AT 


R 


2 


a^ay(^-AT  'S-K 


Xj,(B)  -  X;,(A)  +  ¥^(8)  -  ¥^(4) 
^  2o’a  (B-A) 

X  y 

Y2(B)  -  ¥2(^  -  X2(B)  +  Xj,(^ 
aj  (B-A) 

VarfPkj^l  =  Rj  VarFA]  +  R^  VarfB] 


+  Rj  Var[a^l  +  rJ  VarFa^l 
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